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Abstract. We investigate the evolution of star-forming gas-rich disks, using a 3D chemodynamical model including a dark halo, 
stars, and a two-phase interstellar medium with feedback processes from the stars. We show that galaxy evolution proceeds along 
very different routes depending on whether it is the gas disk or the stellar disk which first becomes unstable, as measured by 
the respective Q-parameters. This in turn depends on the uncertain efficiency of energy dissipation of the cold cloud component 
from which stars form. 

When the cold gas cools efficiently and drives the instability, the galactic disk fragments and forms a number of massive clumps 
of stars and gas. The clumps spiral to the center of the galaxy in a few dynamical times and merge there to form a central bulge 
component in a strong starburst. When the kinetic energy of the cold clouds is dissipated at a lower rate, stars form from the gas 
in a more quiescent mode, and an instability only sets in at later times, when the surface density of the stellar disk has grown 
sufficiently high. The system then forms a stellar bar, which channels gas into the center, evolves, and forms a bulge whose 
stars are the result of a more extended star formation history. 

We investigate the stability of the gas-stellar disks in both regimes, as well as the star formation rates and element enrichment. 
We study the morphology of the evolving disks, calculating spatially resolved colours from ffie distribution of stars in age and 
metallicity, including dust absorption. We then discuss morphological observations such as clumpy structures and chain galaxies 
at high redshift as possible signatures of fragmenting, gas-rich disks. Finally, we investigate abundance ratio distributions as a 
means to distinguish the different scenarios for bulge formation. 
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1. Introduction 



The formation of galaxies is one of the important questions 
in current astrophysical research. With HST and earth-bound 
10m class telescopes it is now possible to study the evolu- 
tion of galaxies with redshift by direct observation. Although 
no complete picture is available at the moment, some facts 
are well established. Surveys undertaken around z ^ 1 
have shown that luminous elliptical and spiral galaxies were 



ISM: structure - Galaxies: abundances - Galaxies: bulges - Galaxies: evolution 



At z ^ 0.5 spiral galaxies and barred galaxies are observed, 
and it seems that by then the full Hubble sequence is in place, 
quite similar to the galaxy distribution today. However, the fre- 
quency of barred galaxies drops sharply beyond z - 0.5 (e.g., 
Abraham et al., 1999). van den Bergh et al. (2002) showed that 
this is not an selection effect. They shifted a local galaxy sam- 
ple to higher redshifts and concluded that most of the barred 
galaxies would still be visible. Thus, the absence of bars at 
higher redshifts seems to be real. 



the present time jBrinchmann et al.l 


19981 iLillv et all Il998 


Abraham et al 1999; 'Abraham', 1 1999 


;lAbraham & Merrifield 


2000; Dickinson. 2000). Beyond z = 


1 .4, there are far fewer 



high-luminosity g alaxies o ff all types compare d to low redshifts 
( lDic kinson''200rf see also Driver et al. 1998). Thus the bulge- 
to-disk Hubble sequence of gala xies appears to have formed at 
redshifts slightly beyond z ~ 1 dAbraham & MerrifieldL Eooot 
iKaiisawa & Yamada 1200 ll) . although part of the stellar popu- 
lations of these galaxies might have formed earlier. 



Already at z ~ 1 around 30% of galaxies are morpho- 
logically peculiar. At higher redshifts many clumpy structures 
and compact objects are observed. Most of these cannot be 
attached to the tra ditional Hubble scheme, but show iiTegular 
morp hologies (e.g.. lAbraham et aLill996tlvan den Bergh et all 
I2OO0I) . The anomalous morphologies observed in these high 
redshift objects cannot be explained through band-shifting ef- 
fects alone, b ecause the iiregula rities persist also in NICMOS 
observations ('Dickinson', l2000l) . probing the visual restframe 
wavelength of these objects. 
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Observations of high redshift galaxies currently provide in- 
formation only about the global properties of these objects. 
Detailed data like stellar metallicity distributions, stellar kine- 
matics, or gas distributions are only available for local galaxies. 
To understand galaxy formation in a consistent picture, mod- 
els must be developed that can be compared with observations 
over the whole observed redshift range. These models should 
be able to explain the properties of distant galaxies as well as 
the detailed data on local galaxies. 

With high resolution cosmological simulations, large 
progress has been made in understanding cosmic struc- 
ture formation (e.g ., Navar ro et al ]_|1996: Mo oreetal.Lll998t 
IJenkins et all l200k .Klvoin et all 1200 L) . However, on galac- 
tic scales these simulations still lack the necessary resolution 
to describe the processes relevant for baryon dissipation and 
star formation. Therefore two approaches for describing galaxy 
formation and evolution have been developed. Semi-analytical 
modeling, based on simple assumptions to describe the bary- 
onic physics and star formation in the dark halos that form 
in the cosmological simulations, has been u sed to analyze the 
global prope rties of g alaxy sa mples (e.g.. iKauffmann et all 
ll993.:.Guiderdoni et al.L199&.Cole et alll2000l) . 

On the other hand, dynamical models, using a subset 
of the cosmological information as initial conditions on 



structure of forming galaxies 


(e 


g..lSteinmetz & Mullet. 


1995, 


Navarro & Steinmet 


7, 1997 


Sommer-Larsen et al. 




1999; 


Williams & Nelson. 


2001,; 


Samland & Gerhard. .2003.). 
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small-scale dynamical models still need to describe star 
formation with a simple parametrization, but they contain a 
much more detailed description of the dynamics, feedback, 
and, in some cases, element enrichment. Thus with these 
models it is possible to predict observable properties through 
the galaxy assembly process, and to compa re directly with 
observations of hi^h redshift galaxies (e.g., IContardo et all 
Il998t IWestera et a l.. '2002: Abadi et al.,'2003'). Because these 
models can be calculated over a Hubble time, one can also 
directly compare detailed present-day characteristics (e.g., 
metallicities and kinematics of stellar populations) with obser- 
vations of local galaxies. Finally, with such models predictions 
can be made for dynamical and stellar population properties 
of high-redshift galaxies, which can be verified by future high 
resolution observations with the next generation telescopes. 

The mass accumulation into dark matter halos is well un- 
derstood in the context of the cosmological simulations and 
can be used a s an input to model g alactic evolution (e.g., 
Ivan den Boschl l2002t IWechsler et all |2D02). Additionally, 
the angular momentum distributions o f the form ing dark 
halos have been c alculated recently ( Bullock et al... .200ll 
IChen & Jingl l2002l) . It is widely accepted, both in the semi- 
analytical approach and the small-scale dynamical models, that 
galactic evolution depends sensitively on the mass and angular 
momentum distribution of the system. 

Much less is known about the processes that govern the 
evolution of the baryons within a dark halo, and how these in- 
fluence the properties of the forming galaxies. Stars in galaxies 
nearby are observed to form in molecular clouds much denser 
than the ambient medium in which these clouds are embedded. 



Most of the kinetic energy of this cloud fluid is in the motions 
of single clouds relative to the bulk flow. This kinetic en ergy 
can be dissipated by inelastic coll isions dL arson, 1969) and 
augmented by supernova feedback ('McKee & Os triker. .19771) . 
Under some conditions the macroscopic cloud system can be 
treated as an isothermal fluid (ICowie et all ll98Cft . Its energy 
dissipation rate is not well-determined, however. One expects 
that it depends on the geometrical structure of the clouds, on 
whether a major part of the dense medium is arranged in fil- 
aments, a nd on their self-gravitating structure and magnetic 
fields (Ki m et alll200lUBalsara et all EoOli). Thus the dynam- 
ics of the macroscopic cloud medium may be more or less dis- 
sipative, depending on the physical conditions, and may well 
vary between galaxies. 

In the present paper we investigate the formation and dy- 
namical evolution of galactic disks, varying the cloud dissipa- 
tion rate. We use the interaction network and two -phase chemo- 
dynamical model of ISamland & GerhardI (l2003[ SG03) to de- 
scribe the assembly of a disk of stars, hot gas, and star-forming 
cold gas. We find that the dynamical stability of the disk de- 
pends sensitively on the cloud dissipation efliciency, here de- 
scribed by the parameter 77coii- For large rjcow, dynamical insta- 
bilities in the gas dominate the evolution, leading to fragmen- 
tation of the disk into a small number of star-forming clumps 
which subsequently merge to form a centrally concentrated 
bulge. For small t^coU, on the other hand, the system forms stars 
until the stellar disk becomes unstable, leading to a stellar bar 
at late times. The diff'erent stability properties in both cases can 
be quantified in terms of the eflFective Toomre Q param e ter for 
the stellar, gaseous, aiid combined system (e.g..lToomre 1 19641 



the stellar, gaseous, aiid combined system (e.g..lToomre 119641 
iTog & Solom'onlll984llmng & Silklll994tlElmegreenlll99.5l) . 



The evolution of a system wit h high dissipation ca lculated 
at higher resolution is described in llmmeh et all (l2003l) . where 
it was shown that morphological and photometrical proper- 
ties of sever al high redshift objects, like the chain galaxies 
(ICowie et al 1 [1995), can be explained by a fragmented disk 
model. 

Several authors have suggested that galactic bulges can 
form by the secular evolution of a galactic disk, driven 
by interstellar gas or stars (e.g.. ICombes & Sandera Il98ll 
IPfenniger & Normanl Il990t iNoguchil Il999l) . Recent observa- 
tions lend some support to these secular evolution scenarios. 
Bulges in late-type spirals show similar properties to their sur- 
rounding disks, in that their light profiles are better fit by 
an exponent i al rath er than an 7?'^"* law (Courteau et all Il996 ; 
Seig ar et all l2002l) . and in their colours dPeletier & Balcells , 
.19961) . 

As shown here, the evolution of star-forming galactic disks 
may take diff'erent routes, depending on whether it is the gas or 
the stars that drives a disk instability. This in turn depends on 
the uncertain efficiency of energy dissipation in the cold cloud 
component. These different routes also lead to different forma- 
tion scenarios for galactic bulges. Thus we suggest here that the 
morphological properties of galaxies may depend not only on 
cosmological variables like different mass or angular momen- 
tum distributions, or infall history, but also on internal physical 
processes during galaxy evolution. 
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In Sect. 12 we describe our model for star-forming disks. 
In Sect. |3] we describe the morphological evolution of these 
disks, depending on the dissipation efficiency of the cold gas. 
Sect. 0] discusses the stability and star formation rates in the 
models, and Sect.|5]describes the properties of the bulges that 
form in the two main evolutionary scenarios. Finally, Sect. |6l 
summarizes our findings. 

2. The Model 

We use a two-phase model for the interstellar medium, consist- 
ing of a hot, low-density phase and a cold cloud medium from 
which stars are formed. The chemical elements and most of the 
energy released from SNell and later SNela are returned to the 
hot phase. However, the cloud vel ocity dispers ion of the cold 
phase is heated by SNell as well jMcKee & Ostriker. 19_77 ). 
We describe this system with a three-dimensional chemody- 
namical evolution code, which combines a hydrodynamical 
grid code for the two phases of the interstellar medium (ISM) 
with a particle mesh code for the stars. See SG03 for more de- 
tails. 

The interactions between the different ISM phases and the 
stars are described in detail in SG03. A difference exists in the 
characterization of the star formation rate (SFR), where we use 
here a simple Schmidt law (>SchmidL.1959.) . 

Psf = Csfp"y (1) 

with a = 1 .5 and Csf c hosen to be consi stent with the star for- 
mation rule derived bv lKennicuttI ( Il998l) . We here adopt a vol- 
ume density star formation threshold pth, converted from the 
surface densit y threshold of 1 - 10 Mopc"^ obtained from ob- 
servations by iKennicuttI lll998h . assuming a disk scale height 
of the order of the spatial resolution of our model. Thus pth - 
1 Mopc^^/SOO pc = 0.002 Mopc"^. This is a lower v alue com- 
pared to other thresholds used in the literature (e.g. iNoguchil 
11999 ^. It ensures that gas does not clump simply because it is 
below a high star formation threshold and cools when there is 
no heating from star formation. 

SG03 did an extensive investigation of the network of pro- 
cesses connecting the different ISM phases in this model. We 
use their values for most of the efficiency parameters that oc- 
cur in this network, as determined through either theory or ob- 
servation. As shown by SG03, this network is strongly self- 
regulating, so that the system is not sensitive to the precise val- 
ues used for these efficiencies. They conclude that the cloud 
dissipation efficiency t^coU is the most uncertain parameter in 
the model. We now describe the relevance of this parameter in 
somewhat more detail. 

The dynamics of the cold cloud medium (CM) is described 
in a statistical way through the hydrodynamic moment equa- 
tions. Energy dissipation and gain are implemented through 
sink and source terms to the moment equations. The internal 
energy describes the kinetic energy of the single clouds rela- 
tive to the bulk motion. The CM cools through reducing this ki- 
netic energy, which in ffiis model is assume d to happen thro ugh 
inelastic collisions between single clouds. iLarso n[ (ll969l) de- 
rived the dissipation rate in a CM consisting of spherically sym- 
metric clouds with constant mass. Using the mass-radius rela- 



tion for clouds given bv lElmegreenI ( Il989l) the internal energy 
change can be written as 

|=-^couP4"'V-^ (2) 

where p is the density of the CM, cr its velocity dispersion, 
and P4 - {Picmlk)l\Q'^ with Picm denoting the pressure of the 
intercloud medium . All constants have been merged to t/coU. 
//coll also contains the deviation of the effective cross section 
for cloud-cloud collisions from the geometrical value due to 
magnetic fields, self gravity, and gravitational focusing. It also 
depends on whether the cold cloud medium is arranged mostly 
in clumpy structures or mostly in filaments, which is not well- 
understood. Thus the precise macroscopic description of cloud 
dissipation on a galactic scale is uncertain, and hence the value 
of ?7coii in eq. |2] It is therefore possible that the dissipation 
efficiency changes from one galaxy to another, and that this 
changes galactic evolution. It is equally possible that interstel- 
lar cloud physics is nearly universal, and that the main driver 
for galactic evolution is the baryonic infall rate per unit area, or 
disk growth time, as .NoguchiCl99 9) has argued. In both cases, 
the local heating-collisional cooling equilibrium of the cloud 
medium is changed, and hence the star formation history. 

To investigate the influence of cloud dissipation on disk 
evolution, we performed a sequence of simulations with dif- 
ferent dissipation efficiencies //coU- The values given in Tabled 
are normahzed to the efficiency t/^jj used in SG03. The range 
of 77coii values corresponds to the maximum variations given in 
Table 1 ofSG03. 



model index 


'7coii ['7f„ii] 


model index 


'7coll 


A 


20 


E 


1/5 


B 


10 


F 


1/10 


C 


5 


G 


1/20 


D 


1 







Table 1. Dissipation efficiencies used in the model sequence, 
normalized to model D, where T/f^y = 0.025 (SG03). 



The dynamical set-up of the models describes the rapid 
formation of a massive galactic disk in a pre-existing ^ stati c 
dark matter halo with NFW-profj l e (.Navarro et al.t 119971) . 
According to ISommer-Larsen et al ] ll2002h the delayed infall 
of the baryonic matter into the relaxed halo can solve the an- 
gular momentum problem arising from ACDM structure for- 
mation simulations. The primordial gas enters the simulation 
volume vertically at |z| = 15.5 kpc, and the infall is uniformly 
distributed over a radius of 17 kpc, with a rotation velocity 
equal to the circular velocity at the infall point. The infall 
rate is 120 Mgyr"' during one Gyr, resulting in a total mass 
of 1.2 ■ 10" M0. The simulation volume has a diameter of 
38 kpc and a vertical height of 31 kpc with a spatial resolution 
around 500 pc. The evolution of a higher-resolution simulation 
of model A is discussed in Im meli et al. (2003) . The higher res- 
olution does not lead to different results, which indicates that 
the outcome of our simulations is not sensitive to the resolution 
used. 
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The chemodynamical model naturally provides ages and 
metallicities of all stars formed over time, as well as ISM den- 
sities and metallicities. This enables us to calculate HST- and 
UBVRIJHKLM-colours of the model, including dust a bsorp- 
tion. T his is done using the method described in Wester a'et alJ 
( I2OO2I) . but we adopt here a three times lower absorption coef- 
ficient. We will use the observable colours to discuss the mor- 
phological evolution. 

3. Morphological Evolution 

3.1. Settling of the Disk 

The infall of the baryonic matter leads to the build-up of a 
gaseous disk. The settling of the disk is shown in Fig.^ which 
shows a cut through the jcz-plane of the density distribution for 
the cold gas component in models B, D and F at different times. 
The energy feedback from the first stars heats up the disk. In 
consequence the settling of the disk takes longer than a free-fall 
timescale, depending strongly on the cloud velocity dispersion 
and hence on the cloud dissipation efficiency rjcow- Since the an- 
gular momentum does not play a role in the settling along the 
rotation axis, the z-component of the velocity dispersion alone 
prevents the clouds from completely falling to the disk plane. 

After 500 Myr the cold gas phase shows a filamentary dis- 
tribution in all three models. This is caused by the local energy 
feedback of the first supernovae of type II (SNell) to the gas. 
This leads to local variations in the velocity dispersion of the 
cloud medium, which act like the analogues of pressure gra- 
dients in the hot gas phase. Thus the cold gas is compressed. 
In addition, in regions with enhanced energy input, the equilib- 
rium between hot and cold gas phases is shifted to decrease the 
density of the cold medium. 

The model with high dissipation (B, first row) settles to a 
thin disk which is almost completely supported by rotation. The 
timescale for the settling is relatively short because the energy 
input from SNell can be dissipated efficiently. The less efficient 
energy dissipation of model D causes the disk to settle more 
slowly. The disk in this model also remains thicker during the 
whole evolution, because the equilibrium between the contin- 
uous energy input from SNell and energy dissipation by cloud 
collisions is shifted to higher internal energies and thus higher 
cloud velocity dispersions. This is even more true in model F, 
where the even less efficient energy dissipation results in a yet 
thicker disk. 

3.2. Face-On Evolution 

Fig.Elshows an overview of the face-on morphological evolu- 
tion of all our models, presented in terms of the V-band rest- 
frame surface brightness calculated from the distribution of 
stars as described in Section 2. Model indices are indicated at 
the left for each row, and time is indicated at the top for each 
column. All images are normalized to the same surface bright- 
ness magnitude interval, which enables us to directly compare 
luminosities of the models. The absolute V-magnitude My of 
the galaxy is also given in each panel. Obviously, the mor- 



phological evolution depends strongly on the dissipation effi- 
ciency. 

At 0.5 Gyr the formation of the disk is under way and 
none of the models shows an asymmetry in surface bright- 
ness yet. Generally, the colder models are brighter in abso- 
lute V-magnitude than their warmer counterparts. At 0.75 Gyr 
the coldest model (A) shows several knots of high density 
with enhanced star formation. In model B a ring like structure 
is visible. 200 Myr later the disks of these two models have 
fragmented, showing several nearby knots. At 1.05 Gyr some 
clumps in models A and B have already merged and fragmen- 
tation is also visible in model C. Whereas at 1.2 Gyr models 
D to G still show a symmetric surface brightness distribution, 
the clumps in models A to C have akeady fallen to the center 
and merged to form a bright bulge. This persists for the remain- 
der of the simulation, together with a clear spiral pattern in the 
disk. Model D develops a bar like structure with spiral arms at 
around 2 Gyr. A very pronounced bar is visible in models E to 
G at around 2.8 Gyr. As can be seen in the last column of Fig.|2l 
these bars become shorter. Bar evolution is discussed further in 
Sect.ES 

The morphological evolution of models A, B, and C is sim- 
ilar in the sense that they all develop a fragmented disk and 
end up with a compact central concentration. Their exists a de- 
pendence on 77coii, in such a way that the more efficiently the 
gas dissipates energy, the more pronounced is the clumping, 
because the lower velocity dispersion in the more dissipative 
models counterbalances the gravitational instabilities less well. 
In all the cold simulations a ring develops. This ring structure 
probably res ults from the form of the rotation curve in these 
models; see ICha & WhitworthI (l2003l) : iGingold & MonaghanI 
( 1983|. It is important to note that the mass enhancement in the 
ring amounts to less than 10% of the total mass in the cloudy 
medium. 

The galaxy models A-C in their fragmented phase resem- 
ble, when viewed e dge-on, the chain galaxies first reported by 
I Cowie et al.1 d 19951 CHS95). These are high redshift galaxies 
observed with HST in the Hawaii Survey Fields, with large 
major- to-minor axis ratios, knotty structures and very blue col- 
ors. CHS95 suggest that chain galaxies in the redshift range 
0.5 - 3 have a mass comparable to that of a present-day galaxy 
and that they represent a new population of galaxies. Our model 
favors the interpretation of O'Neil et al. (2000) that the chain 
galaxies are clumpy disk-like structures seen edge-on. Model 
A, recalculated with higher resolution, is compared to obser- 
vations in detail in Immfili et al. (2003J. The main results are 
that the extraordinary morphology as well as very blue colors 
observed by CHS95 are reproduced well. An investigation of 
the properties of the individual clumps in this model leads to 
typical stellar clump masses of a few 10'' M©. Metallicity dif- 
ferences between individual clumps are not larger than 0.25 dex 
and the influence of individual clumps on the rotation curve can 
by as high as lOOkm/s. 

A similarity in the evolution can also be seen in models E 
to G. In these models the disk remains symmetric for a much 
longer time. At around 2.7 Gyr these models then become bar 
unstable, and a very prominent bar develops in about one dy- 
namical time. Contrary to the fragmented disk models, the ab- 
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Fig. 1. Cut through xz-plane of the cold gas density distribution for models B, D, F. Evolutionary time and model index are 
indicated. 



solute luminosity of the models showing a bar instability does 
not vary much over the simulated time interval. The evolution 
of model D lies somewhere in between these two cases, show- 
ing a bar with spiral arms at around 2 Gyr. 

Comparing models A and B shows that no qualitative dif- 
ferences in the evolution are visible. The same is true for mod- 
els F and G, which indicates that the chosen values of rjcon 
(Table [0 covers the whole relevant range. 

4. Global Properties 

4.1. Stability of the Disk 

As shown in the previous section, our models can be divided 
in two groups according to their dynamical evolution: Models 
A to C show an early fragmentation, whereas in models E to 
G a bar instability occurs at later times. Model D appears to 
be a transition case. In the remainder of this paper, we con- 
fine the discussion to models B, D, F, as representative models 
for the different evolutionary paths. First, we relate these more 
quantitatively to the s tability of the multicomponent disk. The 
Toomre parameter Q (Safranov,'l960'; Toomre,"l964') is an im- 
portant quantity in investigations of disk stability. In a single 
component gaseous disk Q is given through 



(2gas — 



(3) 



where Sgas is the surface density, k the epyciclic frequency, G 
the gravitational constant, and o-g^s is the radial velocity dis- 
persion of the gas. Q = I is a well-defined stability limit for 
such disks; disks with low velocity dispersions and high sur- 
face densities such that Q < I are subject to radial instabilities. 

In a disk consisting of gas and stars the evaluation of Q is 
more difficult. As described above, the interstellar material in 
our model consists of two gas phases. Because the hot inter- 
cloud medium (ICM) only makes up for about 1 % of the total 
gas mass and because it is pressure dominated and therefore 
stable in a Toomre sense, we can neglect the ICM contribution 
to the Toomre parameter. Henceforth, we consider a composite 
Qes for the stellar and the cold gas phases. 



There have been several attempts to describe the stability 
of a two-component syste ms in the context of the Toomre pa- 
rameter .Jog & Solomorn (11984 derived a dispersion relation 
for a two-component system of stars and gas and evaluated 
it numericallv. 'Romeol ll 19921) derived finite thickness correc- 
tio ns for disk systems. A simple approximation was inferred 
bv lWang& 8113(11994 



1 



1 



ega,s Q stai 



(4) 



where Estars is the stellar surfa ce density and cr^,„r ^ is the radial 
velocity dispersion of the stars. lElmegreenl(ll995l) gives a more 
accurate way of calculating the effective Qeff parameter for a 
tw o component disk. He reformulates the dispersion relation 
of I Jog & Solomo nl (11984 in such a way that the derived Qeff 
is completely analogous to the one-phase Q. We use this more 
accurate method as well as the approachof Wang & Silk ( 1994) 
to evaluate Qetf for the two-phase medium. Since the results 
are the same, we will not further distinguish between the two 
methods. 

Fig-Elshows maps of Qgas, Gstar, and two-component Qeff 
for models B and F. The white regions in the maps denote those 
parts of the disk where the number of stellar particles in a grid 
cell in the plane is less than five and no secure stellar veloc- 
ity dispersion and hence no reliable Qstar could be determined. 
These maps show clearly that the Q values can vary substan- 
tially between different parts of the disk, both radially and az- 
imuth ally. 

At around 0.4 Gyr, not many cells reach the stability bound- 
ary just discussed: all Q's have values above unity. Already at 
0.6 Gyr the situation changes in model B. The gaseous compo- 
nent develops regions of instability, and Qeff reveals an unsta- 
ble disk over a large radial range. For comparison, the disk of 
model F is still stable at 1 Gyr in both components as well as in 
the composite QefF- In model B the violent evolution following 
the instability causes a sharp increase in the mean 2eff ■ The in- 
stability converts energy from ordered motions to random mo- 
tions, as well as driving up the SFR and SN heating, which 
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Fig. 2. Morphological evolution of the seven models in terms of V-band restframe surface brightness. The panels show the whole 
simulation area of 38 kpc squared. The model indices are given at the left and time is indicated at the top. Between panels, time 
does not progress in constant intervals; this is to emphasize the interesting evolutionary phases of the models. My is given in 
each panel. 
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indirectly increases the cloud velocity dispersion and hence the 
Toomre parameter. 

At 1.8 Gyr the stellar disk of model F has a central Toomre 
parameter around 1.5. The gaseous phase has a Toomre param- 
eter of at least 2 and therefore is stab le not only to radial but 
also to bar instabilities jPolvachenkol 1 1 9971) . Nevertheless, the 
composite Qeff already approaches values around 1 in the cen- 
ter. Thereafter, geff drops further in the center, causing a bar 
instabiUty at around 2.6 Gyr. The bar is visible in the instabil- 
ity map, but the signature becomes weaker as the bar evolves. 

Globally, both models are stable after their respective in- 
stability, indicating that an equilibrium state has been reached 
and no further large morphological changes will occur (com- 
pare Sect. 15.21 However, model B remains marginally unstable 
in its outer parts, which drives the spiral pattern visible at late 
times (Fig|2l. Consistent with the fact that this model has a 
thinner disk with higher rotational support than model F, the 
outer regions of the disk in model B also have a lower Qeff than 
in model F. 

In summary, the most striking difference between the two 
models is that the instability in model F occurs in the stellar 
disk in the central regions, at comparably late times, whereas 
in model B the instability develops in the cold gas phase over 
the whole of the disk, starting early-on. Thus the different mor- 
phological evolution discussed in the previous section can be 
attributed to the instability emerging in different phases. 

To investigate the time evolution of the Toomre parameter, 
we averaged it over a central circle with radius 2 kpc. The pro- 
gression of this average Q is shown in Fig. |4] Again we see 
that different components trigger the instability in the different 
models. In all models, the instability starts when Qeif - 0.8. 
This value should not be taken as an absolute limit for stability, 
however, since it depends on the radius over which the Toomre 
parameter is averaged. Nevertheless, this shows that 2eff is a 
good tracer of the instability. Also clearly visible is the increase 
of Q after the onset of instability. 

4.2. Star Formation Rate 

The different dynamical evolution along the sequence of mod- 
els has a strong impact on the global star formation history. 
Fig. |5] shows the SFR as a function of time, for models B, 
D, and F. One immmediately recognizes large differences in 
the absolute values as well as in the shape of the SFR. In the 
beginning, the SFR is low because the gas must first fall to 
the Galactic plane and reach the limiting volume density for 
star formation, pth. This takes around 200 Myr. After this the 
SFR increases as more material falls in. The energy feedback 
from the first stars heats the gas and the equilibrium found de- 
pends on the energy dissipation rate of the cold gas phase. The 
less efficiently energy can by dissipated, the higher the cold 
cloud velocity dispersion, and hence the more extended in the 
z-direction the cold gas disk is (Fig. [Q. If, on the other hand, 
the dissipation efficiency is high, the velocity dispersion drops 
and the material settles in a high density disk. Since all models 
have identical mass infall and total mass, the density in the disk 
plane is lower in the thicker disks. Because the SFR follows a 
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Fig. 4. Time evolution of the Toomre parameter for the stel- 
lar component (top), the cold gas component (middle), and the 
composite Toomre parameter (bottom), for models B (black), 
D (dark grey) and F (light grey). The locations where the insta- 
bility emerges are indicated by dashed lines. 

Schmidt law (Eq.Q, the overall SFR will be higher in the mod- 
els with strong gas dissipation. This effect can be seen already 
at around 500 Myr, when the SFR rates of the models differ by 
a factor of two. 

The subsequent shape of the SFR is dictated by the fur- 
ther evolution. In models A to C the cold gas phase of the disk 
becomes unstable, with the stars following the potential pertur- 
bations induced by the gas. The clumping of the gas leads to a 
strong enhancement of the SFR (Fig.|5] top panel). The clumps 
then fall to the center and, due to the efficient angular momen- 
tum redistribution during this phase, a large fraction of the gas 
is located in the center after the clumps have merged. In this 
stage, the SFR reaches a pronounced maximum, which can be 
as high as 200 M^yr"'. 

The warmer models form stars at a much lower rate, after 
the initial rise driven by the infall. Because the gas phase re- 
mains stable, the cold cloud disk from which the stars form has 
no large density enhancements. Star formation then increases 
again when the instability of the stellar disk sets in. The bar 
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time [Gyr] 

Fig. 5. Global SFR (black line) and SFR in the innermost 2 kpc 
(grey line) in models B, D, and F. In all three models the insta- 
bility, indicated by the respective dotted Une, causes an increase 
in the central SFR. 



causes an inflow of gas to the center, causing an increase in the 
central SFR which also refle cts in t he global SFR (Fig.|5l bot- 
tom panel). Indeed, Aauerri! (Il999l) found a strong correlation 
between the global SFR and the barred structure of galaxies. 

As a consequence of their diverse SFR, the integrated mass 
of stars formed in these models also varies strongly (Fig. |6j: 
whereas after 2 Gyr model B has already converted 75% of its 
baryonic mass into stars, the corresponding value for model D 
is around 40% and that for model F is around 25%. 

This illustrates how different the SFR in these models 
are despite their identical mass infall rates. The cooling efli- 
ciency, acting as a process on small scales, thus influences the 
strength and th e shape of the SFR, consistent with results from 
Ivan den BoschI ( l2002l) . It is therefore difficult to constrain the 
mass accretion rate of a galaxy from the star formation history 
alone, without further assumptions. This shows the need for 
self-consistent modelling of galaxies. 

4.3. Metallicity Distribution 

Fig. shows the metallicity distribution of the stars in models 
B and F after time 4 Gyr. Here we divide the stars in the model 
into halo, disk, and bulge components. The halo is defined as 
all stars with a distance to the disk plane of at least 3 kpc. The 
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Fig. 6. Stellar mass fraction of models B (solid line), D (dotted 
Une) and F (dashed line). 

disk component is defined as all stars with z-distance from the 
plane \r,\ < 0.5 kpc and a minimal distance of 2 kpc from the 
center. All stars within 2 kpc from the center are taken to be the 
bulge component. 

There exist clear differences in the metallicity distributions 
of both models, due to their different evolutionary histories. 
However, in both models the metallicity increases from the halo 
to the disk to the bulge. The halo has the lowest metallicities be- 
cause the metal-enriched clouds falls through the halo towards 
the disk plane. When the clouds settle in the disk, the chemical 
enrichment can take place over a longer time-scale. The sub- 
sequent angular momentum redistribution from the instability 
causes the gas to flow inwards. The inflowing gas cannot es- 
cape from the galactic center, and when it is converted to bulge 
stars, it has the highest metallicities. 

Since the settling of the baryons to the plane depends on the 
dissipation rate (Sect. ITU , the enrichment time-scale for stars 
that end up in the halo is longer for model F than for model 
B, and so one expects a more advanced chemical enrichment 
of the halo in model F. Indeed, the halo metallicity distribution 
of the cold model B shows a peak at [Fe/H] ^ -2, whereas 
the peak in the distribution for model F lies approximately one 
dex higher at [Fe/H] ^ -1. The shape of the observed met al- 
licity distribution of the halo of M31 dOurreU et all l200ll) is 
very similar to that of model F. This points to a relatively slow 
star formation process for the M31 halo stars. Bv contrast, the 
metallicity distribution of the Milky Way halo dChiba & BeersL 
I2OOOI) resembles more the distribution created in the cold model 
B, and points to a faster settling of the gas to the disk in the 
early collapse of the Milky Way (see also SG03). 

In the distribution of disk stars, the differences between 
both models are also clearly visible. In model B, the peak at 
relatively low metallicities mirrors the peak in the SFR at early 
times, before SNela play a significant role in the enrichment . 
This model is compared to data from JSdvardsson et alJ lll993h . 
Note, however, that the observed sample contains stars that 
have formed much later than at 4 Gyr where the model cal- 
culation ends. The slower star formation process of model F 
leads to a steady increase of stars with high metallicities; in 
this model the gas is enriched by SNela before a large fraction 
of the stars is formed. 



Immeli et al.: Gas Physics, Disk Fragmentation, and Bulge Formation in Young Galaxies 



11 



Model B Model F 



15 


Halo 


Halo 




E 10 






\ 




m 




^ 5 


Al ■ 


A 









1 1 1 itYin 1 1 rrVffff^%i \ 1 1 iTi i 




15 


■ Disk 1 




- Disk 






10 












z" 


J 










\ 

^ 5 















"L 


iiiTTTrrr^ 1 






1 1 1 1 




h 




Bulge 


Bulge 














\ 
■z: 






J 
















U I . ..I I... I I.... . . I'J I 

-4-3-2-1 ^4 -3 -2 -1 1 



[Fe/H] [Fe/H] 

Fig. 7. Distribution of stars in [Fe/H] for models B (left) and 
F (right), for Halo, Disk and Bulge stars as defined in the text 
(open histograms). Observational data are plotted as the shaded 
histograms (references in the text). 

A similar argument applies to the bulge. Here the high early 
SFR of model B produces many metal poor bulge stars with 
[Fe/H] ^ -0.7. In model F, by contrast, the evolution is slower, 
leading to a steady increase in the stellar metallicity without 
a peak at low metallicities. Tiede & Terndrup ( 1999) reported 
a metallicity distribution for the Galactic bulge with a similar 
shape as produced in model B. On the other hand, results from 
iRamirez e t al. (2000) and Zoccali et al. (2003) point to a more 
metal rich distribution. The formation of the Galactic bulge is 
discussed further in the next section. 

Although clearly the cloud energy dissipation is not the 
only physical process that influences the metallicity distribu- 
tions in galaxies, it does have an important impact on the dy- 
namical, star formation, and chemical evolution of galaxies. 

5. Bulge Formation 

In our series of models we observe two qualitatively different 
paths to bulge formation. In those models where the gas disk 
fragments and develops massive clumps, the bulge is formed 
through t he merging of t hese clumps at relatively early times 
(see also 'Noguchil Il999l) . In the class of models where a bar 
forms through an instability of the stellar disk, the bulge is 



formed from this bar, at comparably late times. This is sim - 
ilar to the scenario propos ed by [Combes & Sander^ ( Il98ll) . 
IPfenniger & Normanl(ll990l) and l^ha et alJ(ll991l) . 

These two routes do not cover the whole range of 
bulge formation processes: further wa ys in which a bulge 
may be formed are th rough mergers dAguerri et all l200lt 
IScannapieco & T isseralEi ppS)) or in the center of an early galac- 
tic collapse (|Eggen et al.[ Il962|) . It is possible that all these 
routes may play a role in the formation of some bulges. Real 
bulges may even form through a combination of these pro- 
cesses, as indeed happens in the self-consistent models of 
SG03. In particular, some nuclear bulges would appear to be 
good candidates for the rapid SFR processes. 

The least well-studied of these bulge formation processes 
is tha t from merging clumps in a fragmenting disk (Nogu chil 
Il999h .l Immeh et al.l (i200 3 ) compare a high-resolution simula- 
tion of model A with observations and conclude that this pro- 
cess may in fact be at work in the chain galaxies found by 
CHS95. Further observations of similar objects will thus be 
highly interesting. 

Fig. |8] shows edge-on projections of the inner regions of 
models B, D, and F in terms of K-band surface brightness. The 
bulges also stand out in the K-band radial surface brightness 
profiles, plotted in Fig.|9] In the following subsections, we give 
some more specific predictions from our models. 




Fig. 8. Edge-on view of the inner 1 1 kpc of models B, D, and F 
(from left to right), in K-band at time 3.8 Gyr. 
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Fig. 9. Radial surface brightness profile in K-band for models 
B (left) and F (right). The violent evolution in model B leads 
to a more massive bulge than the mass transport to the center 
induced by the bar in model F. 

5A. Clump Merging 

As we have seen, in the cold models (A to C) the disk fragments 
in an early evolutionary state and quickly develops several 
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clumps of stars and gas. The clumps then spiral to the center, 
where they merge, causing a massive bulge to form. Because 
this process leads to very high gas densities, the SFR in the 
clumps and, especially, in the subsequently forming bulge is 
very high. This results in a rapid oxygen enrichment of the 
bulge stars, as can be seen from Fig.^| where [O/Fe] is plotted 
against [Fe/H] (top panel for model B). 




Fig. 10. [O/Fe] against [Fe/H] for bulge stars in the three mod- 
els B, D, F. The different transitions from SNell to SNela dom- 
inated abundances are clearly visible. 



In this plot, [O/Fe] at low metallicities corresponds to the 
IMF averaged ab undance ratio of supernovae of type II (SNell) 
( ISamlandill998h . which enrich the ISM on a timescale of or- 
der 10 Myr. The supernovae of type la (SNela), assumed in our 
model to arise from white dwarf binary systems on a timescale 
of around 1 Gyr, influence the chemical evolution of the ISM 
only at later times. The large amount of iron produced by the 
SNela lowers [O/Fe]. The metallicity at which this happens de- 
pends on how far the enrichment has progressed before this 
time, and so the location of the transition in Fig. ^| gives in- 
formation about the star formation time-scale. In the case of 
the fragmenting model B, the short starburst produces a clearly 
defined path through the [O/Fe] -[Fe/H] diagrams, with a bend 
lying around [Fe/H] ^ -0.4. Thus SNell have enriched the 
ISM already to this high metallicity before the feedback from 
SNela plays an important role. 



After the starburst, 70% of the infalling baryonic matter in 
model B has been converted to stars. The fast SF history is thus 
reflected in a large overabundance of a-elements. This is shown 
for the element Mg in Fig. ^2 The distribution of [Mg/Fe] for 
all bulge stars at time 4 Gyr shown in Fig.^Jhas a strong peak 
at [Mg/Fe] = 0.4. There is also a smaller peak at [Mg/Fe] ^ 
-0.2; this value corresponds to the maximum [Mg/Fe] reached 
in a closed box simulation, and is generated by an extended 
low level star formation at late times when the iron enrichment 
from SNela dominates. 

The distribution in Fig.FTTlfor model B looks not unlike that 
for the Galactic bulge stars measured by McWilliam & RichI 
fl994) and shown as the dotted histrogram in Fig. II II Note, 
however, that the observed distribution will depend on the en- 
tire SF history following the simulated early phase. Thus, at 
the moment it is too early to say whether the old Galactic bulge 
could have formed out of an early fragmenting disk, but it is 
clearly a possibility that is worth further investigation. 
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Fig. 11. [Mg/Fe] distribution for bulge stars (r < 2kpc) formed 
until the end of the simulation (at time 4 Gyr), in models B 
(black line) and F (grey line). The dotted histogram shows 
[Mg/Fe] for 1 1 stars measured by McWilliam & Rich ( 1923). 

The evolution of the velocity dispersion of the bulge stars 
is shown in Fig.^J where both the (cylindrical) radial veloc- 
ity dispersion ctr and the vertical component cr, are given. The 
figure shows two main effects: The first, visible in in both com- 
ponents, is a more or less linear rise of the velocity dispersion 
due to the growing mass concentration in the center of the disk. 

The second effect originates in the instability and affects 
only the radial component. In model B, the violent evolution of 
the galactic disk with the final coalescence to a nuclear bulge 
leads to a strong rise in the radial velocity dispersion: when 
the instability sets in, ctr rises sharply from around 30 km/s to 
150 km/s. The corresponding rise in models D and F occurs 
later; see Section l5^ The z-component of the velocity disper- 
sion increases less through the instability. 

At early times, the vertical dispersion in this model exceeds 
(Jr. We believe this is a combination of the infall model used 
and a massive disk building up in the galactic plane. At the 
end of the simulation, the radial velocity dispersion is around 
200 km/ s whereas the z-component lies around 90 km/s. 
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Fig. 13. K-band restframe surface brightness map as mass tracer for the inner 12 kpc of model F. 




Fig. 12. Stellar velocity dispersion of the bulge stars (r < 2 kpc) 
for models B (black), D (dark grey), and F (light grey). The ef- 
fect of the instability on the velocity dispersion is clearly visi- 
ble in the radial component ctr. 



5.2. Bar Evolution 

Fig-Elshows the evolution of model F during the bar instabil- 
ity, through its face-on K-band surface brightness map. Model 
F is representative for the models with lower cloud energy dis- 
sipation, E to G. The bar forms between 2.6 Gyr and 2.8 Gyr, 
i.e., within a time interval of the order of the dynamical time, 
as is a lso observed in N-body simulations dCombes & Sander^ 
The bar instability causes a sudden increase in the ra- 
dial component of the velocity dispersion (Fig. I12> . whereas 
the z-component stays rather unaffected. 

Right after its formation the bar becomes very pronounced, 
and then weakens again. To do a quantitative investigation, we 



made a Fourier decomposition of the mass distribution to derive 
the length of the bar. Fig. [3] bottom panel, shows that the bar 
becomes shorter with time. The bar length evolves from around 
4.5 kpc to around 3.2 kpc within 1 Gyr. The shortening of the 
bar is accompanied by an increasing pattern speed (Fig.[2] top 
panel). The bar developed in the model is a fast bar with its end 
near the corotation radius. The rising pattern speed can be ex- 
plained by the growing mass concentration in the center, since 
all relevant frequencies at a given radius r depend on V-^centrai- 
This evolution is similar in all three models that form a bar. 

Because at the time when the bar forms (Fig. the gas 
content of models E to G is still high, up to 60%, the forma- 
tion of the bar leads to a significant inflow of gas to the cen- 
ter. There the gas can form new stars (compare Sect. 14. 2> . The 
mass accumulation in the center then weaken s the bar, as pre- 
vious l y observed in numeric al simulations (e.g lFriedli & Benzl 
liNorman e tailll996l) . There is also observational evi- 
dence forjheweakening of bars through central concentra- 
tions: IPas et all (l2003l) recently found an anti-correlation be- 
tween central mass concentration and deprojected bar elliptic- 
ity in a sample of 13 barred nearby galaxies. 

As we have already mentioned in the introduction, obser- 
vations in dicate an absence of bars at intermediate and high 
redshifts ( Abraham et all[l999 ). In the models presented here, 
the formation of an extended stellar disk with surface density 
high enough that it can become bar-unstable, takes a compa- 
rably long time. In addition, in these systems with large gas 
fraction, the bar formed is rapidly weakened due to strong gas 
transport to the center Both effects could occur in high red- 
shift disks, as the early disks might have larger cloud random 
motions because of the stronger infall, and would have higher 
gas contents. The two effects thus might explain the absence of 
bars at higher redshifts. 

The metallicity distribution of the stars in the central 
bar/bulge (Fig.llOt reflects the slower enrichment process in the 
hotter models. The bend in [O/Fe] is shifted to lower [Fe/H], 
because of the smaller SFR in these models. The bump visible 
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Fig. 14. Evolution of the pattern speed (top) and bar length 
(bottom) for models E (black), F (dark grey) and G (light grey). 

in this plot for models D and F at [Fe/H] around -0.4 is mainly 
due to the increase of the central SFR caused by the instability. 
The SNell formed in the central regions dominate the feedback 
again and enhance the abundance of a-elements. This effect 
together with the averaging over the central 2 kpc is also the 
source of the scatter in the plots for models D and F. 

We note that, because we use an IMF-averaged yield for 
SNell, we do not find a broad range of [O/Fe] values at low 
metallicities, as would be expected at early times when bulge 
stars are enriched by individual SNell with different progenitor 
masses ( Arg ast et al. . 2000). The scatter seen for models D and 
F in Fig.[lO|comes from the form of the SFR. 

Contrary to the cold model B, the [Mg/Fe] abundance ratio 
distribution of the stars in model F (Fig.llU has only a single 
peak around [Mg/Fe] ^ -0.2. The different maximum metal- 
licities reflect the fact that the SF histories in the models are 
different (compare Fig.|6|l. 

6. Conclusion 

We have investigated a sequence of models for star-forming, 
gas-rich disks. The most important result from these mod- 
els is that galaxy evolution proceeds very differently depend- 
ing on whether it is the gas disk or the stellar disk which 
first becomes unstable, as measured by the respective Q- 
parameters. In our two-phase, star-forming interstellar medium 
description this depends on how efficiently the cold cloud 
medium can dissipate the kinetic energy it gains from dy- 
namical and feedback hea ting. An additio nal important vari- 
able is the infall history llNoguchil [l999l) . which influences 
the gas density and hence the local self-regulation equilibrium 
Jsamland & Ge rhard, 2003, SG03) and the disk stability. 

WTien the cold gas cools efficiently and drives the instabil- 
ity, the galactic disk fragments and forms a number of massive 
clumps. The stellar disk fragments with the gas because of the 



strong gravity of the clumps, which begin to form stars at a 
high rate because of their large density. The clumps then spiral 
to the center of the galaxy in a few dynamical times and merge 
there to form a central bulge component in a strong starburst. 
This scenario is similar to that discussed by Noguch| lll999l) . 
We show that, in this mode of disk evolution, the bulge forms 
rapidly and early; the unstable region of the disk is completely 
disrupted, and the disk must be rebuilt by subsequent infall. 
Because of the starburst origin, many of the bulge stars formed 
in this way have large [a/Fe] abundance ratios. 

On the other hand, if the kinetic energy of the cold clouds 
is dissipated at a lower rate, stars form from the gas in a more 
quiescent mode, while keeping the kinetic temperature high 
enough and the gas density low enough to prevent the gas 
from becoming dynamically unstable. In this case, an instabil- 
ity only sets in at later times, when the surface density of the 
stellar disk has grown sufficiently high. The system then forms 
a stellar bar, which channels gas into the center, evolves, and 
forms a bulge whose stars are the result of a more extended 
star formation history, i.e., have lo wer \afFe]. This scenario re - 
sembles the evo lution described by Combes & Sanders ( 19§j]), 
IPfenniger & Norman C1990.) and Raha et al. (.1991.) . The com- 
parably long formation time for bars and bulges in this mode 
and the weakening of the bar through inward gas flow in the 
still gas -rich disks might expl ain the absence of bars at high 
redshift jAbraham et al .[ ll999h . 

An example of a nearby galaxy in which the fragmenta- 
tion process may be taking place, is the gas-rich, b lue starburst 
galaxy NGC 7673 jHomeier & Gallagheil Il999l) . At higher 
redshifts, a number of morphologically unusual galaxies may 
be in sir nilar evolutionary s tates, such as the so-called chain 
galaxies ilC9wi9 ^t ^l||l^^ ^ and oth er kn otty galaxies seen in 
the HDF hm den Bergh et allll996l) : see llmmeh etaP ( l2003l) 
for further discussion. On the other hand, the abundance ra- 
tios in the centers of nearby bu lges show only a rnodera te 
a overabundance ([a/Fe]^ 0.2; *Pro ctor & Sanson] (^2002*)), 
arguing that these stars were formed from gas enriched by 
more extended star formation histories. The few [Mg/Fe] abun- 
danc e ratio me asurements available for the Galactic bulge 
(iMcWilham & Rich. J394), on the other h and, point towards 
more enhanced [a/Fe]. Consistent with this, |Pii7iaetalJ ll2002l) 
find a similar enhancement in three Galactic bulge fields as in 
bulge globular clusters from integrated light measurements. 

Clearly, mixed star formation histories would also be pos- 
sible, in which a part of the bulge formed in an early starburst, 
preceding a slower build-up of bulge stars from secular evolu- 
tion of the disk. In the multi-phase galaxy evolution model of 
SG03, which follows the formation of a large disk galaxy in a 
growing ACDM dark matter halo, the final bulge is indeed a 
superposition of these two components with, in that model, the 
larger part of the mass in the secular component. Our models 
suggest that the time-scale for secular disk evolution is com- 
parable to or longer than that for Fe enrichment from SNela, 
while the central starburst occurs before SNela become impor- 
tant. Thus further spatially resolved measurements of [a/Fe] 
in bulges, and measurements of [a/Fe] vs Fe and radius for 
Galactic bulge stars will be key for disentangling bulge forma- 
tion history. 
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It is well-known that the properties and evolution of a 
galactic disk depend on global parameters like mass, angular 
momentum, and infall rate. We have shown here that the evo- 
lution may also depend strongly on the physics of the baryonic 
component, in particular, the uncertain energy dissipation rate 
of the cold cloud medium from which most stars are formed. 
This also implies that one cannot simply derive the accretion 
rate onto a galaxy by determining its star formation history. 
Although all the models in our sequence have an identical gas 
infall, we observe very different SF histories from one model 
to the other. In all models, however, the instability causes gas 
transport to the center, and a subsequent increase of the central 
SFR leading to the build-up of a central bulge component. 

We thank the referee, M. Noguchi, for his prompt report, 
the Schweizerischer Nationalfond for financial support of this 
work under grant 20-64856.01, and the Centro Svizzero di 
Calcolo Scientifico (CSCS) for the opportunity of using their 
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